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Abstract 

I In the original holographic dark energy (HDE) model, the dark energy density is proposed to be pde = 3c^M^,L"^, 

with c is a dimensionless constant characterizing the properties of the HDE. In this work, we propose the generalized 
holographic dark energy (GHDE) model by considering the parameter c as a redshift-dependent function c{z). We de- 
rive all the physical quantities of the GHDE model analytically, and fit the c(z) by trying four kinds of parametrizations. 
The cosmological constraints of the c(z) are obtained from the joint analysis of the present SNLS3-i-BAO+CMB+i/o 
data. We find that, compared with the original HDE model, the GHDE models can provide a better fit to the data. 
For example, the GHDE model with JBP-type c(z) can reduce theXmi„ of the HDE model by 2.16. We also find that, 
unlike the original HDE model with a phantom-like behavior in the future, the GHDE models can present many more 
different possibilities, i.e., it allows the GHDE in the future to be either quintessence hke, cosmological constant like, 
or phantom like, depending on the forms of c(z). 
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I. INTRODUCTION 



Since its discovery in 1998 fv], dark energy (DE) 120 has become one of the most popular research 



areas in modern cosmology. Numerous theoretical models have been proposed in the last decade jsMlOll. 
However, the nature of DE still remains a mystery. 

A popular and interesting approach to the nature of DE is considering it as an issue of quantum gravity 
Ini . Since there is no complete theory of quantum gravity now, we usually consider some effective theories, 
in which some fundamental principles are taken into account, such as the commonly believed holographic 



principle [12]. In 1999, based on the effective quantum field theory, Cohen et al. 11311 pointed out that the 
quantum zero-point energy of a system with size L should not exceed the mass of a black hole with the 
same size, i.e., 

L PA < LMj, , (1) 

here L is the ultraviolet (UV) cutoff, which is closely related to the quantum zero-point energy density, and 
Mpi = 1 / ^SttG is the reduced Planck mass. In this way, the UV cutoff of a system is related to its infrared 
(IR) cutoff. It means that the vacuum energy related to this holographic principle may be viewed as DE 
when we take the whole universe into account (its energy density is denoted as p^e hereafter). The largest 
IR cutoff L is chosen by saturating the inequality, so that we get the holographic dark energy (HDE) density 



(2) 



where c is a dimensionless model parameter. If we take L as the size of the current universe, for instance the 
Hubble radius H~^, then the DE density will be close to the observational result. However, Hsu 11411 pointed 
out that this choice yields a wrong equation of state (EOS) for DE. For this reason, Li lisi suggested to 
choose the future event horizon of the universe as the IR cutoff, defined as, 

da' 



Rh = a \ — = a \ 

Jt a Ja 



Ha 



'2 



(3) 



This choice gives a reasonable value for the DE density, and also leads to an accelerated universe. The 
HDE model based on this ansatz has been proved to be a competitive and promising DE candidate. It can 
theoretically explain the coincidence problem lisj . and is proven to be perturbational stable |161 (see |17- 
2 in for more theoretical studies). It is also found that this original HDE model is much more favored by the 



observational dat a 112211 . even compared with other holographic dark energy models with different IR cutoffs 
see Il25l430ll for more researches). 
The parameter c in Eq. ^ plays an essential role in characterizing the dark energy properties in the HDE 
model, e.g. with the value of c being bigger or smaller than 1, the behavior of HDE in the far future would 



be phantom-like or quintessence-like, which means giving very different ultimate fates for the universe (see 
Sec. V for detail). In the previous works, c is always treated as a constant. However, since there are no 
confirmed evidence at present telling us that c shall be a constant, it is worthwhile, and somehow natural, 
to treat this parameter as a variable(a similar idea has been applied on the holographic ricci dark energy in 



13111). motivated by this idea. In this paper, we generalize the original HDE model by treating the parameter 



c as a free function of redshift z, i.e., 

Pde = 3c{zfMliRf. (4) 

Hereafter, we will denote it as "GHDE" (generalized holographic dark energy) model. Not having a theory 
to fix the form of c(z) at present,it is helpful to take some trial functions to characterize it. As a test, we 
consider four parametrizations of c{z) here, i.e., 

1. GHDEl: c(z) = cq + ci^—; (5) 

1 + z 

2. GHDE2 : c(z) = cq + ci ; (6) 

(1 +zr 

3. GHDE3 : c(z) = ^-^ ; (7) 

1 -I- ci ln(l + z) 

4. GHDE4: c(z) = cq + ci - ln2j . (8) 
These choices are inspired by those parametrizations proposed to study the EoS of DE, known as Chevallier- 



Polarski-Linder parametrization (CPL) |32|], Jassal-Bagla-Padmanabhan (JBP) parametrization 113311 . Wet- 
terich parametrization I34I1, and Ma-Zhang parametrization [Ssd, respectively. The original HDE model can 
be recovered when we take ci = in all these four parametrizations. 

In this work, we study the observational constraints for these GHDE models, by fitting the cosmological 



data from the recently released SNLS3 sample of 472 type la supemovae 113611 . the cosmic microwave 
background anisotropy data from the Wilkinson Microwave Anisotropy Probe 7-yr observations 13711 . the 
baryon acoustic oscillation results from the Sloan Digital Sky Survey data release 7 1381, a nd the Hubble 
constant measurement from the Wide Field Camera 3 on the Hubble Space Telescope 13911 . The original 
HDE model is also investigated for a comparison. 

This paper is organized as follows. In Sec. II, we derive the basic equations for the GHDE model. In 
Sec. Ill, by fitting data, we investigate the cosmological constraints on four GHDE models each with one 
form of c(z) listed in Eq. (|5]l. In Sec. IV, we discuss the different fates of the universe in both the HDE and 
GHDE models. At last, we give some concluding remarks in Sec. V. The methodology used in this work is 
listed in Appendix. 
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II. GENERALIZED HOLOGRAPHIC DARK ENERGY MODEL 



In this section, we derive the basic equations for the GHDE model. We will begin with the basic FRW 
cosmology, and then introduce the GHDE model subsequently. 

A. The FRW cosmology 

In a spatially flat isotropic and homogeneous Friedmann-Robertson-Walker (FRW) universe (the as- 
sumption of flatness is motivated by the inflation scenario), the Friedmann equations read 

//2 = p , (9) 



-2Ml,H ^ p + p. (10) 



where H = a/a is the Hubble parameter, and p is the total energy density of the universe, including the 
components of matter p„,, radiation pr and dark energy pde- For simplicity, we will not consider the case 
with interaction between dark matter and DE in this work. So we have energy conservation equation for 
each component, 

p„, + 3Hp,„ = 0, (11) 

Pr + 4Hpr - , (12) 
Pde + 3//(l + W4e)Pde = . (13) 

Combining with Eqs.lfTTT) and (fl^ . The first Friedmann equation Eq. ^ can be rewritten as. 

Ho \ i-a,,(z) 1 ■ ^ ^ 

Here Hq is the Hubble constant, 0.de{z) is the fractional dark energy density, given by, 

^deiz) = — — - . , (15) 

Pciz) ^MjjH^ 

and Q„,o and Q^q are the present values of fractional density for matter and radiation, respectively. From the 
WMAP observations, 0^0 is 

a,o - ^yoii + O.iniNeff) , = 2.469 X IQ-^h^^, Ne/f - 3.04 , (16) 

where Qjo is the present fractional photon density, h is the reduced Hubble parameter, and Ne/f is the 
effective number of neutrino species. 



B. The evolution of GHDE 



Now let us introduce the HDE scenario into the cosmology. Since the original HDE model with constant 
c is a subset of the GHDE model, we can start straightforwardly with the GHDE model. 

First of all, let us derive the effective DE EoS of the GHDE model. Taking derivative for Eq. (|4]) with 
respect to x = Ina, a = -nr, and making use of Eq. Q, we get 



, _dpde ^ ( ,^^dc\ 



Combining Eqs. ([T3] ) and ([TT] ). it follows the DE EoS, 

1 2 ^ ^ 2 {\+z) dc 

Wde ^ -T - T- i^de + T T ■ (18) 

5 5c 5 c az 

Compared with the original HDE model, there is an additional term | ^ due to the redshift dependent of 
c(z). This additional term means that not only the value, but also the differential of c{z) plays an important 
role in the evolution of DE. In the next sections, we will find this additional term make GHDE model 
obviously different with the original HDE model. 

Now let us have a look at the property of the GHDE fraction ^deiz), which influences the expansion 
history of the universe. To do this, we shall solve out the unknown Q.deiz) in Eq. ([T4l ). Directly taking 
derivative for Q.de - c^/{H'^L^), and using Eq. ([3]), we get 

n-2a„(.-,.^.i^), (19) 

where e = -H/H^ = -H'/H. From Eqs. ^ and we have 

3 3 ll \dc\ n^/^ 

e = -il+Wde^de) - o - + -:r]^de - — ■ (20) 

2 2 \2 caxj c 

Thus, with 4- = -(1 + z) J-, we have the equation of motion, a differential equation, for D.de, 



i.a..<,-n„,fi.i:^^-?<ili>$'|. (21) 



dz I + z \ c c dz 

Also, we see an additional term -^^7^^ in the bracket due to the redshift-dependent c in the GHDE model. 
In the data analysis procedure, this useful Eq. (|2TI ) can be solved numerically. And once we have Q.de, the 
Hubble parameter can be obtained directly through Eq. (fT4l ). 



III. OBSERVATIONAL CONSTRAINTS ON THE GHDE MODELS 

In this section, we will place the cosmological constraints on the GHDE models discussed above. Here- 
after we consider four kinds of parametrizitions of c{z) listed in Eq. ^ respectively. 



A. GHDEl: the CPL type 



For this case, c(z) takes the form 

c{z)^co + ci-^ . (22) 
1 +z 

To be more accurate, this parametrization satisfies: z ^ 0, c ^ cq; z — > +oo,c — > cq + ci. i.e. c smoothly 
varied from cq + ci to cq from the past to present. If we assume that the GHDE has positive energy density 
(fide > 0), we can obtain the following constraints on the parameters, 

Co > , CO + Cl > . (23) 

Taking derivatives for Eq. (IHI) . it follows that 

(1 + z) dc Cl 



(24) 



c dz co + (co + ci)z 

Substituting this equation to Eq. ( fTSl ) and Eq. (|2T]) . one can get the EOS and Friedmann equations in the 
GHDEl model. 

In the left panel of Fig. [T] we plot the Icr and 2cr contours for the GHDEl model in the cq - ci plane. It 
is noticeable that the constraints of Eq. (1231 ) (plotted in the black dashed line) is automatically satisfied by 
the results of the numerical simulation. Another feature of this model is that there is a strong degeneracy 
between the parameter cq and ci. The constrained parameter space distributes roughly along the line co + 
Cl - 0. 

In the right panel of Fig. [H we show the evolution of Weffiz) of this model( the w^ff of the original HDE 
model is also showed). Unlike the original HDE model, here we found that the w - -\ line lies in the Icr 
error of the evolution of Weff{z), implying that ACDM model is a subset of this GHDE model. 

B. GHDE2: the JBP type 

The "JBP-typ" GHDE model has the following form of c(z): 

z 

c{z) ^ Co + Cl . (25) 

(1 + zY 

And it follows directly that, 

{\+z)dc ci(l-z) 
c dz co(l + zY + c\z 

From Eq. ( |25l ). we get: z ^ 0, c ^ co; z ^ l,c ^ co + |ci; z ^ +oo,c cq. Like the CPL-type 
parametrization, the JBP parametrization is also unable to describe the behavior of HDE in the future. 




FIG. 1: Left Panel: Marginalized probability contours at Icr and 2cr CL in the cq - ci plane, for the GHDEl model. 
The constraints Eq. (l23b are automatically satisfied. Right Panel: The evolution of Wefj along with z, for the same 
model. The w - - I line and the 1 cr constraint of the original HDE model are also shown. They all lie in the 1 cr region 
of the GHDEl model. 

In this model, the constraint pde > gives, 

Co > , CO + ^ci > . (27) 

In the left panel of Fig. |2l we plot the Icr and 2cr contours for the GHDE2 model in the cq - c\ plane. 
Again, the requirement of pde > is automatically satisfied. The evolution of Weff of this model is plotted 
in the right panel of Fig. |2j Like the GHDEl model, here we also found that both the Weff of the original 
HDE model and the cosmological constant roughly lies in the \cr error of the evolution of w^ff- 



C. GHDE3: the Wetterich type 

The third parametrization considered is the "Wetterich-type" parametrization, having the following form 

of c(z), 

c{z) . (28) 

1 + ci ln(l + z) 

It has the property: z — > 0, c — > cq; z — > +oo, c ^ 0. And it is also straightforward to have 

(1 + z) dc c\ 



c dz 1 + ci ln(l + z) 

The condition p^e > reduces to. 



(29) 



Co > , ci > . (30) 
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FIG. 2: Left Panel: Marginalized probability contours at lo" and 2cr in the cq - ci plane, for the GHDE2 model. 
The constraints Eq. ( |27] ) are automatically satisfied. Right Panel: The evolution of Weff along with z in the GHDE2 
model. The w - -I line and the Icr constraint of the original HDE model are also shown. They are all consistent with 
the GHDE2 model in the Icr CL. 



1.5 



0.5 




-0.6 




FIG. 3: Left Panel: Marginalized probability contours at Icr and 2cr CL in the co - c\ plane, for the GHDE3 model. 
The constraints Eq. dSOl l are slightly violated at small cq. Right Panel: The evolution of w^fj along with z in the 
GHDE3 model. The cosmological constant is slightly deviated at high redshift, while the original HDE model is still 
consistent with this model in the Icr CL. 



The contours of this model are plotted in the left panel of Fig. [3l Here we find that the requirement 
of pde > is slightly violated at the small value region of cq - ci plane. The evolution of Weff of this 
model is plotted in the right panel of Fig.|3] Unlike the GHDEl and the GHDE2 model, here we found that 
both at the very low (z < 0.1) and high (z > 1.2) redshift regions, the behavior of HDE deviates from the 
cosmological constant, while the original HDE model is still consistent with this GHDE model in the Icr 
CL. 
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D. GHDE4: the Ma-Zhang type 



A common short-coming of the above three parametrizations is that they all diverge when the redsliift 
z approaches -1, so we are unable to investigate the future behavior of GHDE in these models. In this 



c(z)^co + ci( ^"P^^^ -ln2| ■ (31) 



section, to overcome this, we consider the following form of c(z), which is proposed in i35h : 

/ln(2 + : 
\ 1+z 

It follows directly that 

z ^ : c ^ Co; z ^ +oo : c ^ cq - ci ln2; z — > -1 : c ^ cq + (1 - ln2)ci , (32) 

^ ^^"'^o + c,(it±^-ln2)- ^ ^ 

We can find that the c{z) of this parametrization is constrained in a finite value even in the future of 
the universe evolution. Thus this model can have a well-defined property in the far future evolution of the 
universe, just as the original HDE model which possesses a constant c. 

In this model, the requirement of pde > is reduced to, 

co>0, co-ciln2>0, cq -i- (1 - ln2)ci > . (34) 

In the left panel of Fig. |4l we plot the Icr and Icr contours for the GHDE4 model in the cq - c\ plane. 
The constraints of Eq. ( [34l ) are automatically satisfied by the numerical simulation results. We also plot 
the evolution of Weff{z) of this model in the right panel of Fig. |4l A most striking property we find in 
this model is that the future behavior of Wgffiz) could be completely different: while in the original model 
the behavior of DE is phantom-like in the future, the GHDE in this model can behave like a quintessence, 
cosmological constant, or phantom. All these possibihties are allowed according to the 2cr constraint of the 
data(see Sec. IV for detailed discussion). 

E. A Summary of the GHDE Models Considered 

In the above subsections, we study four GHDE models with special c(z)s. A brief summary of these 
models are shown in Table HI where the models, model parameters (together with their best-fit values and 
Icr uncertainties), a.ndx'^ ■ s are given. To make a comparison, the results of the original HDE model and the 



ACDM model (see Ref. 14011 ) are also listed. Here some nuisance parameters, such as h, a and /3 mentioned 



in the appendix, are not shown since they are not model parameters with significant meanings. 
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FIG. 4: Left Panel: Marginalized probability contours at Icr and 2cr CL in the cq - ci plane, for the GHDE4 model. 
The physical conditions, Eq. ( l34b . are automatically satisfied. Right Panel. The evolution of Weff along with z in the 
GHDE4 model. The cosmological constant and the Icr constraint on the original HDE model all lie in the Icr region. 
Diff'erent from the previous three GHDE models, the future evolution of Weff is also plotted. 



TABLE I: Fitting results for the GHDE models. 



Model 




Co 


C\ 


^ mm 


ACD 








424.911 


HDE 


261+0021 

'^■^"^ -0.019 


65+*"" 





424.855 


GHDEl 




^•■^^-2.86 


294+3.16 


423.369 


GHDE2 




'^•^-'-0.19 


65+°-24 

'-'•"-'-0.93 


422.696 


GHDE3 


262+0 026 


1-18-0.63 


50+'°^ 

"•-'"-0.62 


423.735 


GHDE4 




97+0-31 

'-'•^'-0.45 


1 TO +0.47 


423.407 



From the Xmin^ we can see that all these models can provide a nice fit to the observational data. Especially, 
the GHDE2 model showed that the inclusion of the extra parameter can reduce theXm^ of the original HDE 
model by 2.16. 

By fitting the observational data, we also get the constraints for those four special c(z) forms, which are 
plotted in Fig. |5] It is found that both the best-fit and the constraint regions of c(z) are different among the 
four GHDE models. As mentioned in Sec. IIIBI plays an essential role in characterizing the evolution 
of DE in the GHDE models. Thus the difference can be interpreted by the different forms of ^^s in these 
models. 

We also find that, although all the c(z)'s best-fits of these four GHDE models distinctly deviate from 
c's Icr region of the original HDE model, the Icr constraint of the original HDE model lies well in the Icr 
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FIG. 5: The evolution of cs along with z for the four considered GHDE models. For having a loosely closed probabiUty 
contours in the cq - c\ plane(see Fig.[T]and Fig.|2]i, thus the upper bounds of c{z) are missed in both the GHDEl and 
GHDE3 models. 



region of each GHDE models considered in this work, so we conclude that there are not enough information 
to determine the form of c{z) in the GHDE scenario based on the current cosmological observation data. 

IV. THE FATE OF THE UNIVERSE IN THE GHDE SCENARIO 

The fate of the universe is always a fascinating issue[41]. In this section, we discuss the fate of the 
universe in the GHDE scenario. First we study the fate in the original HDE model as a comparison. 

A. The fate in the original HDE model 

As studied in Sec. JIJ the behavior of HDE is essentially determined by the parameter c. In the original 
HDE model, the parameter c is proposed as a constant, thus the Eq. (fTSl ) reduces to, 

1 2 I 

Wde = -^- T' y^de ■ (35) 
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FIG. 6: Left Panel: Marginalized probability contours at Icr and 2cr CL in the f2,„ - cq plane, for the HDE model. We 
obtain c < 1 in the Icr CL, implying that the future behavior of HDE is phantom like. Right Panel: The evolution of 
Weff along with z for the HDE model. The cosmological constant with w = -1 is also plotted. In the far future when 
z approaches -1, we have Wejf < -1 in the 2cr CL. 



From the discussion in Sec. Ill A[ we can find that in the far future, the DE will be the dominated composite 
of the universe, i.e., z ^ -1, Q.de 1 . From Eq. (l35l) . we have, 

1 2 

Wde ^ - :r' z ^ -1 . (36) 

3 3c 

From Eq. ( [36b the relation between c and the fate of the universe can be seen clearly. If c = 1 , in the far 
future HDE will behave Uke a cosmological constant; if c > 1, then always we have w > -1 and HDE 
behaves as a quintessence; and if c < 1, we will have w < -1 in the future, leading to a phantom universe 
with big rip as its ultimate fate. 

Based on the joint analysis from the SNLS3+BAO+CMB-i-//o data, in the left panel of Fig. |6j the Icr 
and 2cr contours for the HDE model in the Q.,„ - cq plane are plotted. In the 2cr CL, we obtain 

n,„ = 0.265!°:EJ36, c-0.65!°:I^, (3?) 

showing that c < 1 is significantly favored by the data, which implies that, for the orginal HDE model, 
there will be Wgff < -1 in the far future (see the right panel of Fig.|6ll. Thus, in the orginal HDE model the 
behavior of HDE is phantom like at last, and the universe will probably end up with a big rip i4lll . 



B. The fate in the GHDE models 



In the GHDE scenario, the parameter c is not a constant but a function of the redshit z- This makes 
the variety of the HDE behavior increases, and thus let the GHDE model present many more different 
possibilities for the fate of the universe. 
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We can focus on the GHDE4 model for a typical demonstration. From Fig. |4l it is interesting to see that 
various possibilities emerge even in Icr region: the behavior of GHDE in the future evolution can be either 
quintessence-like or phantom-like, which means the big rip can either happen or not. 

Inspired by this fact, one may wonder what forms of c(z) can make the HDE processing a purely 
quintessence-like, phantom-like, or even cosmological-constant-like behavior. This issue is to some ex- 
tent equivalent with another issue, that is, how to reconstruct a given form of DE component in the GHDE 
scenario. 

For a most general form of DE component, we can describe it by the DE density function f(z), 

Pdeiz) = Pdeofiz) . (38) 

Combining with Eqs. ([3]) and we have 

*) = v^wu) = r ^ (39) 

1-i-z J_i £(z') 

So indeed we can utilizing the GHDE model to reconstruct a DE component with arbitrary from, once given 
its energy density function f{z). 

As an example, let us consider the case that the DE component with a constant w, i.e. f{z) = (1 -i-z)^'^^''""'^ 
Neglecting the small Q.^, thus 

dz' 



:(z) = V(l-n^o)(l+z)i+3" 



(40) 



Vn^oCl + Z')' + (1 - i^m0)(l + Z')'(^+") 

A demonstration is shown in Fig. |7J We take the value Q,„o = 0.25 and w = -0.8, -1.2,-1, which 
correspond to a purely quintessence-like, phantom-like, and cosmological constant behavior, respectively. 
By numerically integrating the Eq. ( |40l ). we get the evolution of c(z) for these ws. Especially, the case 
w = - I gives the well-known ACDM model. Thus, by choosing different forms of c(z)s, we can reconstruct 
a DE component with various behaviors. 

As already discussed in the above section, that we can not confirm the form of c(z) based on the cun^ent 
cosmological observation data. So, in the GHDE scenario, the fate of the universe remains a mystery. 



V. CONCLUDING REMARKS 



In this work, we investigate the generalized holographic dark energy models by considering a redshift- 
dependent c{z). In all, we consider four parametrizations of c(z), including the CPL type, JBP type, Wet- 
terich type and Ma-Zhang type parametrizations. By performing the joint data analysis, it is shown that all 
these models can give a nice fit to the data. Among them, the JBP-type GHDE model can even lead to a 
reduction in the . of 2. 1 6. 
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FIG. 7: The evolution of c along with z that corresponding to different value of constant w. especially w = -1 is the 
cosmological constant case. 

From the fitting results we also find that the generalization of c{z.) increase the variety of the HDE model. 
However, the evolutions of Wgffiz) in these models are all consistent with that of the original HDE model in 
lo" CL. That means, we do not see clear evidence for the evolution of c(z) based on the current observational 
data. 

The fate of the universe is always a fascinating issue. Unlike the original HDE model predicting a 
probably big-rip in the far future, we find that the GHDE model, with the form of c(z) undetermined by the 
current data, lay this issue back to a mystery. Thus for the GHDE model, big-rip or not, it continues to be a 
question. 

At the end, we hope that the method of considering the parameter c of the HDE as a redshift-dependent 
function may provide us some clues to understand the DE problem. In all, more studies about this scenario 
are needed in the future. 

Acknowledgements 

We would like to thank Xin Zhang for helpful discussions. This work was supported by the NSFC under 
Grant Nos. 10535060, 10975172 and 10821504, and by the 973 program (Grant No. 2007CB815401) of 
the Ministry of Science and Technology of China. 








14 



TABLE II: Relationship between number of free model parameters rip and Ax^. 



Tip 






1 


1 


4 


2 


2.30 


6.18 


3 


3.53 


8.02 


4 


4.72 


9.72 


5 


5.89 


11.31 


6 


7.04 


12.85 


7 


8.18 


14.34 


8 


9.30 


15.79 


9 


10.42 


17.21 



VI. APPENDIX: OBSERVATIONAL DATA AND METHODOLOGY 

In this work, we adopt the statistic to estimate the model parameters. For a physical quantity ^ with 
experimentally measured value ^obs, standard deviation cr^ and theoretically predicted value ^th, the;^^ takes 
the form 

= -2 ■ (41) 



The total is the sum of all;k'|s, i.e. 



= 24 • (42) 



One can determine the best-fit model parameters by minimizing the total x^- Moreover, by calculating 
Ax^ = x^ ~^mm' determine the Icr and the 2cr confidence level (CL) ranges of a specific model. 

Statistically, for models with different np (denoting the number of free model parameters), the Icr and 2cr 
CL correspond to different A^^. In Table HH we hst the relationship between Up and Ay^ from ?ip = 1 to 
Hp = 9. 

In this work, we determine the best-fit parameters and the Icr and 2cr CL ranges by using the Monte 
Carlo Markov chain (MCMC) technique. We modify the publicly available CosmoMC package 114211 and 
generate 0(10^) samples for each set of results presented in this paper. We also verify the reliability and 
accuracy of the code by using the Mathematica program I^bI- 



For data, we use the SNLS3 SNIa sampl e Bo] , the CMB anisotropy data from the WMAP7 observations 
3?! . the BAO results from the SDSS DR7 |38], and the Hubble constant measurement from the WFC3 on 
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the HST 13911 ■ In the following, we briefly describe how these data are included into the^^ analysis. 



A. The SNIa data 



Here we use the SNLS3 SNIa dataset released in [36]. This combined sample consists of 472 SN at 



0.01 < z < 1.4, including 242 SN over 0.08 < z < 1.06 from SNLS 3-yr observations 13^1, 123 SN at low 
redshifts |^| 



|4J], 93 SN at intermediate redshifts from the SDSS-II SN search 14^, and 14 SN at z > 0.8 
from HST |47l]. The systematic uncertainties of the SNIa data were nicely handled Js^. The total data of 
the SNLS 3 sample can be downloaded from Ii48i1 . 
Thex^ of the SNIa data is 



(43) 



where C is a 472 x 472 covariance matrix capturing the statistic and systematic uncertainties of the SNIa 
sample, and Ani = iAb - n^mod is a vector of model residuals of the SNIa sample. Here niB is the rest- 
frame peak B band magnitude of the SNIa, and mmod is the predicted magnitude of the SNIa given by the 
cosmological model and two other quantities (stretch and color) describing the light-curve of the particular 
SNIa. The model magnitude mmod is given by 



5 logio ^^iCzhei, Zcmb) - ais - \) + /3C + M . 
Here Dl is the Hubble-constant free luminosity distance, which takes the form 

-Zcmb 



2)L(Zhel,Zzcmb) = (1 + Zhel) 



Jo 



(44) 



(45) 



E{z') ' 

where Zcmb and Zhei are the CMB frame and heliocentric redshifts of the SN, s is the stretch measure for 
the SN, and C is the color measure for the SN. a and j6 are nuisance parameters which characterize the 



stretch-luminosity and color-luminosity relationships, respectively. Following 13611 . we treat a and fi as free 
parameters and let them run freely. 

The quantity M in Eq. ((44)) is a nuisance parameter representing some combination of the absolute 
magnitude of a fiducial SNIa and the Hubble constant. In this work, we marginalize M following the 



complicated formula in the Appendix C of 13611 . This procedure includes the host-galaxy information 14911 
in the cosmological fits by splitting the samples into two parts and allowing the absolute magnitude to be 
diff'erent between these two parts. 

The total covariance matrix C in Eq. (I44l ) captures both the statistical and systematic uncertainties of 
the SNIa data. One can decompose it as Ii36ll . 



C — Dcfat + Cstat + C 



sys 



(46) 
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where Dstat is the purely diagonal part of the statistical uncertainties, Cstat is the off-diagonal part of the 
statistical uncertainties, and Csys is the part capturing the systematic uncertainties. It should be mentioned 
that, for different a and /?, these covariance matrices are also different. Therefore, in practice one has to 
reconstruct the covariance matrix C for the corresponding values of a and yS, and calculate its inversion. For 
simplicity, we do not describe these covariance matrices one by one. One can refer to the original paper 



Bol and the public code i48h for more details about the explicit forms of the covariance matrices and the 
details of the calculation ofx^^^- 



B. The CMB data 



Here we use the "WMAP distance priors" given by the 7-yr WMAP observations 13711 . The distance 
priors include the "acoustic scale" Za, the "shift parameter" R, and the redshift of the decoupling epoch of 
photons z*- The acoustic scale Ia, which represents the CMB multipole corresponding to the location of the 



acoustic peak, is defined as ll37h 



Ia = (1 +z*) 



rsiz*) 



(47) 



Here Z)a(z) is the proper angular diameter distance, given by 



(48) 



and rsiz) is the comoving sound horizon size, given by 

-l/(l+z) 



V3 Jo 



da 



a^H{a) sjl + (3Qi/4ay)a ' 



(49) 



where Q.b and Qj are the present baryon and photon density parameters, respectively. In this paper, we 
adopt the best-fit values, % = 0.02253/2"^ and Qy ^ 2.469 x I0~^h'^ (for Tanb ^ 2.725 K), given by the 



7-yr WMAP observations 1371] . The fitting function of z* was proposed by Hu and Sugiyama 15011 



(50) 



where 



gi 



0.0783(afo/i^)-"-^^^ 
1 + 39.5(afo/j2)0-763 ' 



g2 



0.560 



1 -F21.1(Qfo/j2)i-8i ■ 



(51) 



In addition, the shift parameter R is defined as 15111 



R{z*) = ^n„,//2(i + 2JDa(z*) . 



(52) 
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This parameter has been widely used to constrain various cosmological models 15211 . 
As shown in Is?], the ^2 Qy^g ^j^ta is 



(53) 



where Xj = {Ia,R,z*) is a vector, and (C(^L„),-, is the inverse covariance matrix. The 7-yr WMAP observa- 



tions i37h had given the maximum likehhood values: Za(z*) = 302.09, R{z*) = 1.725, and z* = 1091.3. The 



inverse covariance matrix was also given in 13711 



2.305 29.698 -1.333 
29.698 6825.27 -113.180 
-1.333 -113.180 3.414 



(54) 



C. The BAO data 



Here we use the distance measures from the SDSS DR7 13811 . One effective distance measure is the 



Dviz), which can be obtained from the spherical average 15311 

Dviz) ^ 



(l+zfDliz) ^ 



H{z) 



1/3 



(55) 



where Da(z) is the proper angular diameter distance. In this work we use two quantities cIqi = 
rs{Zd)/Dv{0.2) and Jo ss = rs(Zd)/^y(0.35). The expression of is given in Eq.(|49l). and Zd denotes the 
redshift of the drag epoch, whose fitting formula is proposed by Eisenstein and Hu IsJ] 

1291(Q,„/j2)*'-25i 



Zd 



1 + 0.659(Q^/l2)0.828 



where 



= 0.313(a™/i2)-"-4i'^[l +0.607 (O^/j^)"-^''^], 



Following 



where 



b2 - 0.238{nmh^f^^\ 



we write the for the BAO data as. 



Api = pf'^ - Pi, pf'^ = d^^ = 0.1905, ^ = d^^^ = 0.1097, 



and the inverse covariance matrix takes the form 



30124 -17227 
-17227 86977 



(56) 

(57) 
(58) 

(59) 

(60) 

(61) 
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D. The Hubble constant data 



The precise measurements of Hq will be helpful to break the degeneracy between it and the DE param- 



eters i55h . When combined with the CMB measurement, it can lead to precise mesure of the DE EOS w 



1560. Recently, using the WFC3 on the HST, Riess et al. obtained an accurate determination of the Hubble 



constant 13911 



Hq = 73.8 ± 2.4km/s/Mpc, 



(62) 



corresponding to a 3.3% uncertainty. So the^^ of the Hubble constant data is 

\2 



h - 0.738 
0.024 



(63) 



E. The total 

Since the SNIa, CMB, BAO and Ho are effectively independent measurements, we can combine them 
by simply adding together the;^f^ functions, i.e.. 



2 _ 2 2 2 2 

Xau - Xsn '^Xcmb '^Xbao ^Xh- 



(64) 
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